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Objectives 


After working through this unit you should be able to: 


(i) locate the local maxima and minima of a function of one variable; 

(11) determine the greatest or least value of the images of a function 
of one variable; 

(111) determine the behaviour of a function of one variable near the 
zeros of its derived function; 

(iv) sketch a rough graph of a given (simple) function; 

(v) express the equations of simple surfaces using Cartesian co-ordinates 
in three dimensions; 

(vi) represent a function of two variables as a_ three-dimensional 
surface ; 

(vii) define the partial derivatives of a function of two variables, and 
evaluate these derivatives in simple cases; 

(viii) determine the equation of a tangent plane to a given surface at a 
given point; 

(ix) determine the local maxima and minima of a function of two 
variables; 

(x) distinguish between local maxima and minima of a function of 
two variables, x and y, by determining the behaviour of curves of 
intersection between the surface and planes through the point P, 
where the tangent plane to the surface at P is parallel to the xy- 
plane. 


Note 


Before working through this correspondence text, make sure you have 
read the general introduction to the mathematics course in the Study 
Guide, as this explains the philosophy underlying the whole course. You 
should also be familiar with the section which explains how a text is 
constructed and the meanings attached to the stars and other symbols 
in the margin, as this will help you to find your way through the text. 
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Glossary 


Terms which are defined in this glossary are printed in CAPITALS. 


CONTOUR LINE 


FUNCTION OF 7 
(REAL) VARIABLES 


LOCAL MAXIMUM 
(MINIMUM) OF A 
FUNCTION OF ONE 
VARIABLE 


LOCAL MAXIMUM 
(MINIMUM) OF A 
FUNCTION OF TWO 
VARIABLES 


OVERALL MAXIMUM 


(MINIMUM) 


PARTIAL DERIVATIVE 


PARTIAL DERIVED 
FUNCTION 


SADDLE POINT 


STATIONARY POINT 
OF A FUNCTION OF 
ONE VARIABLE 


v1 


A CONTOUR LINE is a line on a SURFACE marking a 
particular level; for example, a subset of points on 
the surface with equation z = F(x, y), with co- 
ordinates (x, y, k), where k is a constant. 


A FUNCTION OF n (REAL) VARIABLES is a function 
whose domain is a subset of 


Rx kK x =- x &, 
en] 


ntimes 


and whose codomain is R. 


If f is a FUNCTION OF ONE VARIABLE, with domain 
A, and there is a positive ¢ such that f(x) < f(c) 
(f(x) > f(c)), where ce A, for all 
xEAnl[c—eé,c + ¢€], then f has a LOCAL MAXIMUM 
(MINIMUM) at c. 


If F is a FUNCTION OF TWO VARIABLES, with domain 
A, and there is a positive ¢ such that F(x, y) < F(a, b) 
(F(x, y) > F(a, b)), where (a, b) € A, for all 


(x, YE AM {(x, y):(x — a)? + (y — bY < 87}, 
then F has a LOCAL MAXIMUM (MINIMUM) at (a, 5). 


The OVERALL MAXIMUM (MINIMUM) of a (real) func- 
tion is the greatest (least) value of the set of images 
of all the elements in the domain of the function. 


If F is a FUNCTION OF TWO VARIABLES, x and y, then 
the PARTIAL DERIVATIVE OF F WITH RESPECT TO 
THE FIRST VARIABLE (1.€. x) at (x, y) 1S 
; S Fis.> Ay} — FO) 
(x, y) = bm —————-, 
h~0 h 
if this limit exists. 
The PARTIAL DERIVATIVE OF F WITH RESPECT TO 
THE SECOND VARIABLE (i.e. y) at (x, y) is 
. F(x, y + k) — F(x, y) 
1%.7)* 
k~+0 k | 


if this limit exists. 


If F is a FUNCTION OF TWO VARIABLES, x and y, then 
the PARTIAL DERIVED FUNCTION OF F WITH RESPECT 
TO THE FIRST VARIABLE (i.e. x) is the function: 


F’, :(x, y)—> Fi, y), 


and the PARTIAL DERIVED FUNCTION OF F WITH 
RESPECT TO THE SECOND VARIABLE (i.e. y) is the 
function: 


F’, :(x, y)+ —> F(x, y). 
A SADDLE POINT 1S a STATIONARY POINT of a FUNC- 


TION OF TWO VARIABLES which is neither a LOCAL 
MAXIMUM nor a LOCAL MINIMUM. 


A STATIONARY POINT OF A FUNCTION f OF ONE 
VARIABLE, With domain A, is a point c € A such that 


ite} =@. 
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STATIONARY POINT 
OF A FUNCTION OF 
TWO VARIABLES 


SURFACE 


TANGENT PLANE 


A. STATIONARY POINT OF A FUNCTION F OF TWO 
VARIABLES, with domain A, is a point (a, b)e€ A 
such that F(a, b) = F,(a, 6) = 0. 


A SURFACE is the geometric figure consisting of 
those points whose co-ordinates, x, y, z, satisfy any 
equation such as | 


z = f(x, y) 
or 
Fix 9. 2)-= 2 


Alternatively, x, y and z may be given by three 
parametric equations of the form: 


x eo, a), y= Heese). 
z= AG. o} 


where G, H and K are FUNCTIONS OF TWO VARIABLES, 
r and « (called parameters). 


The TANGENT PLANE to a SURFACE at a point P 
(on the surface) is the plane which is such that each 
line in the plane which passes through P is a tangent 
to the surface at P. 


Vil 
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Notation 


The symbols are presented in the order in which they appear in the text. 


:& The derived function of the function f. 
— The set of positive real numbers. 
In The natural logarithm function. 
SS The derived function of the function f’, that is, the second 
derived function of f. 
e e = exp(1) = 2.71828... 
F A function of two (real) variables is usually denoted by a 
capital letter. 
F(x, y) The image of (x, y) under the function F. 
Ee The Cartesian product set, R x R. 
R" The Cartesian product set of R™” and R,R x Rx... x R. 
ee 
vn terms 
(x, y, Z) The ordered triple, having x as its first element, y as its 
second element, and z as its third element. | 
F; The partial derived function of the function 
(x, y)-—> F(x, y) (x, y)e R®) 
with respect to the first variable, x. 
F(a, b) The partial derivative of F (given above) with respect to the 
first variable, x, at (a, b). 
F; The partial derived function of the function 
(x, y)—> F(x, yy) (ye R’) 
with respect to the second variable, y. 
F‘,(a, b) The partial derivative of F (given above) with respect to the 
second variable, y, at (a, b). 
OF 
= Alternative notation for F(x, y) (see above). 
S(a, b, €) Fhe set: 
{(x, y):( — a)? + (y — bY? < 27}. 
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15.0 INTRODUCTION 


Many problems in both pure and applied mathematics are concerned 
with maximum or minimum properties of some sort. For example, at 
what angle should a missile be fired in order to give the maximum range? 
What is the largest area which can be surrounded by a given length of 
fencing? What is the shortest path between two points on a given surface? 
Problems of this kind are sometimes called optimization problems, and 
some of them can be attacked systematically using calculus. 


In this unit we are mainly concerned with the problem of determining 
the greatest and least values attained by the images of a given function, 
and the elements in the domain to which these images correspond. Such 
problems occur quite frequently in practical situations, but, as one would 
expect, the functions which arise in realistic cases tend to be rather formid- 
able. For this reason the “‘practical’”’ situations which we discuss here are 
grossly simplified, and, as in our first example, we sometimes make no 
pretence of realism at all. However, the techniques which we develop 
are often used to discuss non-trivial systems, and it is the wide variety of 
applications which makes these methods so useful. 


We have already seen some methods for solving problems of this kind in 
Unit 6, Inequalities, but the discussion there applied mainly to linear 
functions: here the functions are more general. 


Essentially, the aim of this unit is to bring to your attention some tech- 
niques and ideas which can be developed into rigorous methods for the 
study of optimization, and to point out those areas which will need a 
deeper treatment in later years. 


We begin by discussing functions of one real variable, by which we mean 
real functions (whose domains and codomains are R or subsets of R). 


The second part of the text is devoted to a similar treatment for real- 
valued functions of two real variables (functions which map R x R, or 
subsets thereof, to R). Before tackling this we discuss a little three-dimen- 
sional co-ordinate geometry which we need to know, and which, in any 
case, is of interest in its own right. 


You will probably have most difficulty with sections 15.2 and 15.3, 
mainly because these sections deal with three-dimensional situations. 
For this reason the television programme is devoted almost entirely to 
the geometric notions which we shall need in three dimensions: the 
ideas of planes, surfaces and so on. 


13.1 OPTIMIZATION OF FUNCTIONS OF ONE 
VARIABLE 


15.1.1 A Cautionary Tale 


A mathematically minded Chancellor of the Exchequer once had what he 
thought was a wonderful idea for a new sort of tax. He would levy a tax 
on the amount of overtime that people worked. The idea seemed to be 
superb from every point of view. For their own good he would be deterring 
people from working long hours, and it was merely incidental that a 
large sum would be raised annually for the Exchequer. Thus his humani- 
tarian principles and the requirements of his office would be satisfied 
together. 


This is what he intended to do. For the first hour worked in excess of 
forty hours, the worker would be required to pay 5p, for the next hour 
10p, and so on. What could be easier? 
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A fellow called Fred earned 50p an hour in a factory, and he was puzzled. 
How many hours should he work in order to earn the most money? 


For 40 hours work he would get £20.00; 
for 41 hours work he would get £20.45; 
for 42 hours work he would get £20.85; 


and so on. 


After calculating the total pay for 


all the possible numbers of work- Hours worked Total pay (£) 


ing hours from 40 to 55, he could 40 50.00 
see that it was a waste of time to 

work for more than 49 hours per = tas 

ack 42 20.85 

43 21.20 

44 21.50 

45 Zi.7? 

46 21.95 

47 — 22.10 

48 22.20 

49 2225 

50 2225 

5] 22.20 

52 22.10 

53 21.95 

54 21.75 

55 21.50 


His wife looked over his shoulder (as wives will) as he worked all this out, 
and promptly told him he was a fool. Obviously he should not work for 
more than 49 hours, because the rate at which he earned money went 
down from 50p an hour in the 40th hour by steps of 5p an hour until in 
the 50th hour the rate was down to zero, after which his rate of earning 
was negative and he would simply be giving money away. 


Total Pay 

a ee ee See eee ee ee ee ee eee ee ee ee eee Se ee ee SS ee ee ee 
ZIi@iGiloiei#2#isizf2ifezisti£goil €£f ai é fis 
a ee ee ee eee! et ee Se ee, ee ee ee ee ee ee ee Se ee Se Se ee ee ee ee 
oo; oj;jo;]oj;jsewo;]@ese;]soe;@ea;o;o; oo; oo}; 09; @ 
2 SS ee ee eae eS ee SS ee Ste SS fa Fe =e Se 
GSiatiaqaiarha&iaiaQaiaiate@aite@eiai; @ias; a) s¢6 
o;jw li ocotiws]s o;ywmit] o& wr orrwmo;r ao; wm ;oi;]w 
Si SiS (Bi Sin tate +S Se ae 

e 

ae * % 
& H e 

® H z 

' 
¢ ' e 

' 

' 

e H 

t 

f 

& H 

H 

i] 

e H 

H 

' 

t 

1 

! 

i 

' 

! 

{ 

i] 

1 

' 

' 

t 

' 

H 
Hours Worked 


40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 


Everything proceeded satisfactorily for a time with Fred working a 
regular 49 hour week, but after a while the Chancellor felt that it was time 
to add a few refinements to keep the tax inspectors on their toes, and to 
keep his name before the great voting public. 


What if someone were to work for some fraction of an hour? The Chan- 
cellor had been keenly following the Open University’s Mathematics 
Foundation Course, and was anxious to try out his new-found knowledge 
of finite differences. It was quite a simple matter for him to formulate 
the rule for calculation of ‘“‘Overtime Tax”: 


‘The tax deducted shall be 3(x — 40)(x — 39) new pence, where x is the 
total number of hours (or fraction thereof) worked in one week.” 


He didn’t find it difficult to verify that the formula gave him the same 
result as before when x was any whole number greater than 40. 


Fred then had the problem of finding the value of x which gave him the 
maximum pay. He was fairly certain that it was between 49 and 50, but 
to determine it exactly was simply beyond him. All he could say was that 
his pay would be: 


50x — 3(x — 40)(x — 39) new pence per week. 


Total Pay 


Graph of 
xr2>50x - 3(x-40)(x-39) 


x 
40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 


Luckily Fred never had to solve the problem, for the very next day the 
Chancellor caught a cold and could not go to work for a week. At the end 
of that time he was presented with a bill for £39.00 from the local tax 
office. He promptly changed the law. 


This story brings out some of the important points of the first section 
of this text, in which we wish to optimize functions of one variable. 
Before the introduction of the Chancellor’s amendment, it might seem 
to be a sound practical technique to calculate the total pay for each 
possible number of hours worked, and then choose the largest value. Fred 
chose to do this from 40 to 55 hours, but why stop there? He could be 
sure that his pay wouldn’t be more than £22.25 after 73 hours, say, but 
what would he have done if the tax law had been more complicated? 
The method of looking at the image of every element in the domain of the 
function under consideration (in this case the function which maps hours 
worked to total pay) will only work when the domain is finite. Even then, 
it may involve a considerable amount of calculation if the domain has a 
large number of elements. 


When the domain of the function is not a finite set, such as the interval 
of real numbers in our example, then we must think of some other tech- 
nique. The answer lies in a comment from Fred’s wife : ““Stop working once 
your rate of earning 1s zero, because then you have reached your maximum 
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wage.” There is a little more to the technique we shall develop, but 
essentially that is the idea. As another example on the same line of thought, 
you reached your maximum height once your rate of growth became zero, 
in other words, when you stopped growing. 

In this first section, then, we shall use the concept of rate of change to 
solve optimization problems for functions of one variable. Clearly we 
are going to need some results concerning rate of change from Unit /2, 
Differentiation I. 


Exercise 1 


A supermarket sells chocolate biscuits at 20p a packet and sells 1200 
packets weekly. The manager estimates that for every penny by which 
the price is reduced, sales will increase by 300 packets a week. The lowest 
possible price is 12p, which is what the supermarket pays the wholesaler. 
At what price should the manager sell the biscuits in order to maximize 
his profit? = 


15.1.2 Using the Derivative 

To introduce a method of optimization which uses the derivative, let us 
first look at a fairly simple function. 

Example 1 


What are the greatest and least values of the images of the function f, 
where 


fixe x* —FTx* + Bx + 1 (x € [0, 5])? 


First we look at the graph of f : 


f(x) Graph of 


x R> x3—7x24+8x+ 


Domain 


Imagine the tangent line AB moving along the curve from the point where 
x = 0 to the point where x = 5. The slope of this line is initially positive, 
becomes negative, and is positive again when we reach x = 5. At two 
intermediate points the line is parallel to the x-axis (it has zero slope) 
and the graph shows that these are the points at which f(x) takes its greatest 
and least values in the interval [0,5]. Remember that f’(x) is the slope of 
the tangent at x. If 


f(x) = x° — 7x? + 8x + 1 (x € [0, 5]), 
then we know from Unit 12, Differentiation I that the slope at x is given by 
f'(x) = 3x? — 14x + 8. 


The values of x for which f’(x) = 0 are the two solutions of the quadratic 
equation: 


3x? — 14x + 8 = 0, 


gi Bae Bre ee 


Exercise 1 
(S minutes) 
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Example 1 


(See RB14) 


(See RBS) 
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which are 


and x = 4. 


BOTS 


=, Gee 
The greatest and least values of f(x) in the interval [0, 5] are therefore 
f(@) = 344 and f(4) = -15 


respectively. 


* f(x) 


Graph of 
xX ee x3-7x2+8x+1 
Tangent horizontal 


Domain 


tui FLU ee 


Tangent horizontal = 

Exercise 1 Exercise 1 
: : (4 minutes) 
(1) Find the greatest value of g(x), where 
2 eee oe (x € [ — 2, 2]). 
(11) Find the greatest rectangular area which can be enclosed by a fence of 
length 100 metres. x 
This technique has to be used with a certain amount of care on some Discussion 
occasions. The following examples illustrate the sort of difficulties which 
can occur. 
Example 2 Example 2 
What is the greatest value of 
Fiee=x — Tx? + 8x41 in [0, 6]? 
You may well say that the answer is 334, as in Example 1. But, if so, how 
do you explain the fact that f(6) = 13? 7 
The apparent contradiction is explained when we examine the graph of 
the function, which shows that the greatest value of f(x) in [0, 6] occurs 
when x = 6. 
f(x) 
Graph of 
X Re x3-7x2+8x+1 
Domain 
= (continued on page 7) 
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Solution 15.1.1.1 Solution 15.1.1.1 


Suppose that the supermarket manager fixed a price of x new pence; 
his profit on each packet would then be (x — 12) new pence. He would 
have reduced his sales price by (20 — x) new pence, and therefore 
increased his sales by 300 times (20 — x) packets per week. The total 
number of packets sold per week would be 


1200 + 300 (20 — x), 
so that his profit would be (1200 + 300(20 — x))(x — 12) new pence 
per week. — 


In order to determine the selling price which gives the greatest profit, we 
need only construct the following table: 


ee a __ —————————— 


Selling price, x (p) Profit per week (£) 

12 0 

13 33 

14 60 

5 81 

16 96 

17 105 

18 optimum 108 greatest 
selling profit 
price 

19 105 

20 96 


EE ———e ee — 


There are other ways of solving this problem, as we shall see later. 


Solution 1 Solution 1 


(i) Since x? > 0 for all x, the greatest image is g(0) = 4. To illustrate the 
method we are introducing, we also sketch the graph of g: 


g(x) 


The slope of the tangent to the curve at the point x 1s 
g(x) = —2x 


This slope is zero when x = 0. From the graph of g, we see that 
(0) = 4 is the greatest of the set of images of f. 


FM 15.1.2 


(11) > 


50 — x 


Let x be the length in metres of one side of the rectangle. The area of 
the rectangle is 


x(50 — x) m7, 
SO we can express the area of the rectangle by the function: 
g:x'— x(50 — x) (x € [0, 50}). 
Then 
2'(x) = 50 — 2x 
and 
g(x) = 0 when x = 25. 


With this value for x, the rectangle is a square, and the required area 
is 625 m?. g 


(continued from page 5) 


Example 3 Example 3 
What is the least value of f(x), where 
7 (xe R)? 


This function (which you first met in Unit /, Functions) is called the 
modulus function; it has the following graph: 


Graph of 
xX > |x| 


x 


The difficulty in this case is that f is not differentiable for all values of x. 
We saw in section 12.2.1 of Unit 12, Differentiation I that 


f(x) = +1 ifx > 0, 
and 


f(x) = -1 ux < ©; 


but if x = 0, then the limit: 
TEE bw + h) — f(x) 


h~>0 h 


does not exist, and so f is not differentiable at x = 0. 


Unfortunately in this case it is this very point which is of most interest, 
for we can see from the graph that it is at this point that f(x) takes its 
least value. = 


The following strategy will enable us to deal with most problems (and 
certainly all the problems of this type which occur in the Foundation 
Course) without having to draw graphs. (It is nearly always useful, 
however, to draw a rough sketch of the graph when this is not too 
difficult.) There are cases in which things can still go wrong, but they 
will not concern us here. 


A Strategy for Finding the Greatest and Least Values of the Images of a 
Well-behaved* Function 


Given the problem of finding 
the greatest and least values 


of f(x) where f is continuous and 
has domain [a,b] 


Find the values of x other Find the values of x 
than a and b for which f‘(x) for which f’(x) = 0 


does not exist 


Evaluate f(a) and f(b) 


Evaluate f(x) at these values of x 


Take the least of all 

these values of f(x) 
This is the 

LEAST VALUE 


Take the largest of all 
these values of f(x) 


This is the 
GREATEST VALUE 


Sometimes it is impossible to specify the greatest value of the images of a 
function, simply because there is no such value. Consider the function: 


te eee (x € R). 


In this case, there are always elements in the domain for which f(x) takes 
values greater than any fixed number you care to name. (Notice that we 
may not say ‘“‘the greatest value is infinity” because infinity is not a number, 
and by the words “greatest value” we imply that we are looking for a 
number.) 


* By a “well-behaved” function we mean a continuous function which is differentiable at all 
but a finite number of points. in its domain. It is possible to find a function which is con- 
tinuous at all points in its domain, but not differentiable at any point in its domain, and 
to which the above strategy does not apply. (For the definition of a continuous function, 
see Unit 7, Sequences and Limits I.) 
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Exercise 2 Exercise 2 
(5S minutes) 


Find the greatest and least values (whenever possible) of the images of 
the following functions, and sketch the graph of each function. (See RB14) 
(i) h:x+—300(36x — 288 — x’) (x € [12, 24}). 
(Compare this with Exercise 15.1.1.1.) 


= 1 
(li) f:x-—x + - (xe R*) 
< 
= 16 
2 aan a ice R*) 
x 


Solution 2 


(1) You will find that the supermarket manager’s profit function (Exercise 
15.1.1.1) can be rearranged to give the same rule as the function which 
we are discussing here: 


h -x-——> 306 (36x — 288 — x*) (xe[12, 24)). 
(The domains of the two functions are different.) We have 

h'(x) = 300 (36 — 2x) 
so that h'(x) = 0 when x = 18. This gives the greatest value of h(x), 
which is 10 800p. The least value of h(x) in the given domain is 0, 


and this value occurs at both endpoints of the domain, that is, at 
x = 12 ant x — 724. 


h(x) 


x > 300 (36x -288- x2) 


10 000 


5 000 


= 1 
(i) f(x) =x + 


so that 


1 
PS 


and therefore 


{f@®)=0 4 x= Se 


However, only the value x = 1 is in the domain of f: On this occa- 
sion it helps to sketch a graph of f, and we can easily do this if we 
notice the following: 


= 
(a) x + —is never zero for xe R* ; 
x 
Ze 
(b) x + = is very nearly equal to x when x is large, because the value 
= 
of — is then very small; 
x 
= 1 
(c) x +— is very nearly equal to — when x is small, because then 
x x 


= : | 
— is large and predominates over x. 
x 


Using the above facts we can draw the following sketch: 
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The least value of f(x) is f(1) = 2; there is no greatest value. 
1 
(ii) efx) =x? + = 


so that 
16 
g'(x) = 2x — oe. E 
and therefore 
2'(x) = 6 if 2° — 16-= 0, that is, if x = 2. 
Once again we can draw a sketch if we notice that: 


(a) x? + —isa little greater than x* when x is very large; 


1 
(b) x? + — is a little greater than — when x is very small (but not 
x 
zero); 


(c) g(x) > 0 when xeER™. 


Using the above facts, we can sketch the following graph: 


The least value of g(x) is g(2) = 12; there is no greatest value. 
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15.1.3. Local Maxima and Minima 


Although we often wish to find the greatest (or least) value of the images 
of a function, there are occasions when it is useful to be able to find points 
(like x = 4 in Example 15.1.2.2) where there is a sort of “minor peak’”’ 
on the “‘side of the mountain’’; this is after all a first step towards the 
greatest value. In order to be more precise, we formulate the following 
definitions. 


‘ Let f be a given real function with domain 4A; let ce A. 


If there is a positive number ¢ such that 

f(x) < f(c) for all xeAnlc—ec+é] 
then we say that f has a local maximum at c*. 
If there is a positive number ¢ such that 

f(x) = f(c) forall xeAn[c—ect+el 


then we say that f has a local minimum at c. 


To distinguish these values from the “‘greatest” and “‘least”’ values that 
we have been discussing, we shall call the greatest (or least) value taken 
by the images of a function its overall maximum (or minimum). 


How do the above definitions apply in the context of Example 15.1.2.2? 
The function: 


fix-—> x? — 7x* + 8x +1 (x € [0, 6]) 


has a local maximum at 4. If we take ¢ = 4, say, then f(x) < f(4) for all 
x € (5, 1]. 


f(x) 


Graph of 
f: x R> x3- 7x2+8xt1 


This function also has a local maximum at 6, because for all x in A = [0, 6] 
close to 6 we have f(x) < f(6). For example, taking ¢ = 0.1, the set 
An [c—6,c + €] becomes [5.9, 6], and f(x) < f(6) in this interval. 


Graph of 
f: xe> x3- 7x2+8x+1 


200 001 sa ta et aD oy eae oe co TNA 
es cs ee a hhc ce hn 


* We shall also say, for example, “‘(c, f(c)) (or simply “‘c’’) is a local maximum’”’, when it 
would be more correct to say “‘(c, f(c)) is a local maximum point”’. 
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15.1.3 


Definitions 


xk *® 


Definition 1 


Definition 2 


Definition 3 


xx 


Discussion 


x * 


In this case we know that f(6) is in fact the overall maximum value of 
f(x) for x € A. The reason for taking Am [c — é,c + é] in our definition 
is that we are only interested in that part of the interval [c — ¢,c + é] 
which lies in A. If c is not an end-point of the domain, we are interested in 
the behaviour of f close to c on both sides of c: but, for example, if c is at 
the right-hand end of the domain, we need only look to the left of c. 


Effectively, we are saying that a point which is a local maximum is an 
overall maximum in its immediate surroundings, and similarly for a 
local minimum. Speaking very roughly, if it rains on the graph of a 
function, the puddles collect around the local minima, and the water 
runs away to the overall minimum when the puddles overflow. A local 
maximum would be a suitable place to stand in a flood, but an overall 
maximum would be preferable. 


Of course, an overall minimum is also a local minimum: but a local 
minimum need not be an overall minimum. The following graphs show 
some of the various possibilities. Each of them corresponds to some 
function f. 


f (x) 


Overall maximum 


Overall minimum 


Domain 


f(x) 


Overall maximum 


Overall minimum 
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f (x) 


Overall 
maximum 


Local 
maximum 


Local 


ry 
Overall minimum 
minimum 


Domain 


Exercise ] 


Mark the local and overall maxima and minima on the following graphs: 


(1) 


Domain 


(11) 


1 
Domain 
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Exercise 1 
(5 minutes) 


(iii) 


a a 


Domain 


We call points x such that f’(x) = 0 stationary points of f. A stationary 
point is thus simply a point on the x-axis where the tangent at the cor- 
responding point on the graph is parallel to the x-axis. 


If we wish to locate local maxima (or minima) of a function, it would 
seem a sound idea to first locate the stationary points. However, there are 
unfortunately two complications. 


We have already seen that a local maximum (or minimum) of a function 
can occur at a point which is not a stationary point (in other words where 
the slope of the graph is not zero) either because the function is not 
differentiable at that point, and “‘slope” is meaningless, or because the 
point occurs at an end-point of the domain. We shall overcome this 
difficulty by considering only functions which are differentiable at all 
points of their domains, and by restricting our search for the local maxima 
and minima of such a function to points which are not end-points of the 
domain, and then examining the end-points as a separate issue. 


There is a second complication which is more serious: a stationary point 
may be neither a local maximum nor a local minimum, as we shall see in 
the next example. 
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Definition 4 


Discussion 


x* 


(continued on page 17) 


Solution 1 


(i) 


(11) 


(iii) 


f(x) 


Local 
minimum 


Overall 
maximum 


Overall 
minimum 


ee ee ee a ee ee ne ee ee 
Seiad) 


f(x) 


Local 
minimum 


Domain x 


Overall maximum 


IJU.»UCCns Overall minimum 


i 
Domain 


Local maximum 


Local minimum 


Overall 
minimum 


FM 15.1.3 


Solution 1 


FM 15.1.3 
(continued from page 15) 
Example 1 Example 1 
Consider the function 
7X (x € R) 


which has the graph: 


Graph of 
xX ee x3 


We know from Unit 12, Differentiation I that 


fs 5 ie Bi 


so that 

f'0) = 0. 
and therefore f has a stationary point at 0. But f has neither a local maxi- 
mum nor a local minimum at 0. & 


How then can we distinguish between local maxima, local minima and 
stationary points which are neither? 


15.1.4 Two Useful Methods 


In this section we shall describe two methods for determining the nature 
of stationary points. 


Method One 


Let f be a real differentiable function: 
7a (x € A). 


If f’(x) <0 for all x e S, where S © A, then we say that f(x) is decreasing 
in S. 
If f(x) > 0 for all x e S, where S © A, then we say that f(x) is increasing 
in S. 
If there is a point c € R such that 

Jie <0 4 X= 
and 


foj > -# *Se 


then f(x) is decreasing in {x:x < c} and increasing in {x:x > c}. Clearly, 
f(c) 1s the overall minimum of f(x). 


f(x) 


* iy 
C P 
%, é 
> ®@ 
é 5 
c y 
%, » 
o 


oe 


If there is a point ce R and a positive number ¢ such that 
foo <0 8 c-—t< 2 Se 

and 
Fixt>O 1 Geax cere 


then, although we cannot say anything about overall maxima and 
minima, we can be sure that f(c) is a local minimum of f(x). 
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Technique 


zx=** 


Definition 1 


Definition 2 


Exercise 1 


For each of the following functions determine the subset of the domain 
in which the images of the function are increasing and the subset in which 
they are decreasing. 


fix x’ (x e[—5, 5]) 
(ii) 6:x-— ax” — xx + 1n3 (x € R) 
3 
(iii) g:1 5 - 6 Se # 
Exercise 2 
If 
3 


t 
+ 7 - tei 2 33 


find a positive value of ¢ for which 


(i) g(t) is increasing in ]|1, 1 + e]; 
(ii) g(t) is decreasing in [1 — e, 1[; 
(iii) g(t) > g(1) for te(1 —¢, 1+ ¢]. = 
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Exercise 1 
(3 minutes) 


Exercise 2 
(3 minutes) 
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Solution 1 


Solution 1 
(i) f(x) = 2x, 
which is negative when x < 0 
and positive when x > 0. 
' 


f (x) 


= Domain 


The required subsets are: 
[—5, O[ for f decreasing; 
]0, 5] for f increasing. 
(11) @(x) = 2nx — f, 
which is negative when x < 4 
and positive when x > 4. 


The required subsets are: 
{x:x <4} for decreasing; 
{x:x > 5} for @ increasing. 


(ili) g(t) = 1? — 1, 
which is negative when t? < 1 (i.e. when —1 <t <1) 
and positive when t? > 1 (ie. when t < —1 or t > 1). 


g(t) 


s 
(7) 
& 
®@ 
Cy 
(¢) 
=, 
~ 
S 


“eM 
é, 
eo, 
” 
8 
ce) 
S. 
i) 
Lo} 
@ 
NS 
& 
a 
2 
fe) we 
e) 


OT EERE DES Ss ORR THCERNI WF OT REY TE OY AM OAH ARN DM MOPAR FAIRE EA BSOUYHI on, 


Tk 
| NOI leetertrtestesterietietetetoen 


The required subsets are: 
]—1, 1[ for g decreasing; 
[—2, —1[ u ]1, 3] for g increasing. e 


Solution 2 


For (i), any positive ¢ less than or equal to 2 will do. 
For (ii), any positive ¢ less than 2 will do. 
For (iii), any positive ¢ less than or equal to 2 will do. = 


Method Two 


You may be familiar with another technique which involves the second 
derived function of the function under consideration. This technique does, 
however, require a certain amount of care, not so much in its application, 
which is often very straightforward, but in the conclusions which you draw. 
Unfortunately, most students entirely ignore the method which we 
have just discussed once they are introduced to this second technique. 


We know that the derived function, f’, can be used to study the rate at 
which the images, f(x), of a function f are changing. Similarly, f” can be 
used to study the rate of change of f’(x). In our previous method it was the 
fact that f’(x) changed sign at a local maximum or minimum which was 
important, and it is a short step to the question: “Can we predict a change 
in the sign of f’(x), if we are given f”?”’ 

To find an answer to this question, suppose that we are given a real 
continuous differentiable function : 


gx 2} (x € R), 


and that g’ is also a real continuous function. (We shall say in a moment 
how g is related to f.) Suppose further that g(c) = 0 and g'(c) > 0 for some 
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Main Text 
kkk 
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ce R; then, since g’ is continuous, there must be an interval [c — ¢,c + ¢] 
in which g’(x) > 0. It follows that g(x) is increasing for x E[c — 6c + ¢€]. 
(As shown in the figure below, g(x) may be increasing outside this interval 
as well, but that is immaterial to the argument.) 


ee te ce eee aes ce os ad a te a in en ms A ed le 


{ 
: 
4 
i] ’ 
i f 
g (x) is increasing 
in this interval 


We see from the graph that 

g(x) < g(c) i c-2=< x <<¢G 
and 

g(x) > g(c) fof <x est 


(Don’t forget that one of our assumptions is that g(c) = 0.) 


The useful piece of information which we are seeking follows from the 
above if we now assume that g = f”’. 


Notice that we require g(c) = 0, and this implies that 
Fey = 9%, 


in other words, c is a stationary point of f. We also require that g’(c) > 0; 
since g’ = f”, this implies that 


fc > 0. 


With these conditions we were able to conclude that g(x) < g(c) if 
c—ée<x<c, and g(x) > g(c) if c< x <c +e. Translated to give a 
result for f, this becomes 


f(x) < 0, a = = oe 
and 
f(x) > 0, ifec<x<ct+e. 


This means that there is a local minimum at c. 


If originally we had taken f"(c) < 0, then our final conclusion would be 
that there is a local maximum at c. 


(We assumed for convenience that the domain of f was R, but the same 
results are true for any function which has an interval as its domain.) 
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Classification of Stationary Points Using the Second Derivative 


Local 
maximum 


What should be entered in the blank circle? 
A Few Words of Warning 


Negative 


f'(c) = 0) 


Find the stationary points 
of f (values c such that 


Positive 


What can we say if f‘(c) = Oand f"(c) = 0? It is very tempting to say that 
f has neither a local maximum nor a local minimum at c, but this is wrong. 
._The following examples should make the point clear. 


Example 1 


Consider the following three functions, each with domain R: 


fiz 


f(x) 
f(x) = x* 
f(x) = 4x° 
f'(0) = 0 
f(x) = 12x? 
f"(0) = 0 


From the graph, 
f has a local 
minimum at 0. 


gx | —.x* 


g(x) 
g(x) = 1 — x* 
g'(x) = — 4x? 
(0) =0 
gi = =12 
g°(0) = 0 


From the graph, 
g has a local 
maximum at 0. 


h:x-— x3 


h(x) 
h(x) = x? 
Ory = 3<* 
h'(0) = 0 
"(x)= 62 
h'(0) = 0 


From the graph, 

h has neither a 
local minimum nor 
a local maximum 
at 0. 
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Discussion 
xx 


Example | 


What can we enter in the blank space on our flow chart corresponding 
to the case when f’(c) = 0 and f"(c) = 0? The answer is that the space 
is best left blank, for, given only this information, we can say nothing 
except that the tangent to the graph at c is horizontal. To specify the nature 
of the stationary point, we require more information. There are more 
powerful tests for classifying stationary points using higher derivatives, 
but we leave these until later. (You may find it interesting to try to con- 
struct such a test for yourself.) Remember that Method One works 
even when f"(c) = 0. 


Although calculus is a wonderful tool, it isn’t a substitute for common 
sense. A little concentrated thought will occasionally go a long way, 
as you can see in the following example.* 


Example 2 
Find the overall minimum value of 
g(x) = ((x* + 2) + x*(3 — x*))? (x € R). 


If your first thought is: ““Differentiate, and to the devil with the subtleties”, 
then we admire your single-mindedness, but not your common sense. 


The following solution is much simpler. Simplifying, we get 


a(x) = (x* + 2 + 3x* — x*p 


= (2 + 3x7)’. 
Since x* > 0,2 + 3x? takes its least value, 2, when x = 0: hence the overall 
minimum value of g(x) 1s 4. z 
Exercise 3 


Find the stationary points of the following functions, and classify each 
of them as a local maximum, a local minimum or neither. 


(i) f:x-—> x? — 6x* +9x +6 # (xeR,), 
ny 2x (xe R*). 


* See also Exercise 15.1.2.1, part (i). 
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Discussion 


x * 


Example 2 


Exercise 3 
(3 minutes) 


15.1.5 Examples and Exercises 


You may omit this sub-section at a first reading if you wish, since it 1s 
intended only to consolidate the ideas introduced in the first section of 
this text. On the other hand, you may feel that you need more practice, 
in which case you should work through a number of examples and exer- 
cises until you feel confident enough to continue to the second and third 
sections. 


Be careful not to spend too much time on any particular example or 
exercise; if there is a point which you don’t understand, it is often a good 
idea to read on a little further. 


Example 1 


Of all the cylindrical tin cans (with lids) which enclose a volume of 1000 
cm?, which requires the least metal? s 


Solution of Example 1 


Let S be the surface area of the can in cm7; then, if the metal is of uniform 
thickness, the amount of metal required is proportional to S. Let h be the 
height, and r be the radius of an end-face measured in cm. (We have 
chosen the variables h and r, but the reason for this choice is not obvious. 
Choosing suitable notation is an essential prerequisite to mathematical 
solution, and we refer you to Polya* (page 134) for a discussion of the 
matter.) 


Since the volume enclosed is 1000 cm*, we have 
mr*h = 1000. 
Also, S is given by 
S = 2nr* + 2nrh 
Our problem is to choose values of h and r satisfying Equation (1) which 


will minimize S. 


Equation (2) could be used to define a function (r, h)+—— S, with domain 
R* x R*, but in our problem r and h are not independent, and so we 
cannot choose (r, h) at random from R* x R”. By eliminating h, we can 
obtain a mapping 


jf ise—45 rer) 


*G. Polya, How to Solve It, Open University ed. (Doubleday Anchor Books 1970). This 
book is the set book for the Mathematics Foundation Course; it is referred to in the text 
as Polya. 
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Further Examples 


Example I 


Equation (1) 


Equation (2) 


(continued on page 26) 


Solution 15.1.4.3 


(i) f'(x}-= 3x*.—12x + 9 
= 3(x — 3)(x — 1). 
Thus f'(x) = 0 when x = 1 and when x = 3. 
f(x) = 6x — 12 
and 
f°) = = 6, 
which is less than 0, giving us a local maximum at x = 1. 
f"(3) sia 6, 
which is greater than 0, giving us a local minimum at x = 3. 
= | 
- (ii) Hix} = & x : +Inx 
=1+Inx 
; 1 
Thus h’(x) = 0 when Inx = —1, that is, when x = = 
| 
h"(x) = -, 
x 


— a = 
which is greater than 0 when x = -, giving us a local minimum at 
e 
this point; it is, in fact, an overall minimum. = 


(continued from page 25) 


and then use the methods which we have developed to find the value of r 
which minimizes S$ = f(r). Once r is determined, the corresponding 
value of h follows from Equation (1), and hence we know the dimensions 
of the required can. 


1000 es 
Solving Equation (1) for h, we obtain h = — then substituting for h 
r 


in Equation (2) gives 


S = 2ar* + 2nr|—> 
mr 


ze 


1000 
= 2 pe + “ 


This equation defines a function: 
1000 
- oe (re R*). 
mr 


free 


(Notice that we cannot have r = 0, for both practical and mathematical 
reasons.) 
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Differentiating f, we get 


1000 
f(r) = 2x] 2r — “4 
2nr? — 1000 

= 2 —— 


so f(r) = 0 when 2zr* = 1000. 
The only stationary point of f 1s at 
10 


f= 
J/2n 
Since 
10 
f(r) < 0 when 0 < r < —— 
J/2n 
and 


10 
<2 
it follows from Method One (page 18)* that f(r) takes its overall minimum 
value when 


f(r) > 0 when ‘> 


10 

Vm 

Substituting this value for r into Equation (1), we find that 
20 

San 

so that the tin which requires least metal has h = 2r. 

The graph of f looks like this: 


r= 


Domain 


* You may like to show this by calculating f”(r) and then substituting r = , if you feel 


Pn 


that Method Two is easier. 
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Example 2 Example 2 
Sketch the graph of the function 
P xt {x — PG + 2) (x € R). 


Solution of Example 2 


\ 


(i) Obviously P(x) = 0 
ix = tory = —Z 
and P(x) = 8if 
x = 0. So these 
points are easily 
plotted on the graph. 


(11) When x is a very large 
positive number, 
P(x) is very large and 
positive. When x is 
a very large negative 
number, P(x) is very 
large and negative. 
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(iii) Differentiating P, we 
find (after a little P(x) | 
manipulation which 
we leave as an 
exercise) that P(x) = 
(x — 1)(x + 2)*(5x 4+ 1), 
so that there are 
stationary points of 
P at 1, —2 and — 4. 

It seems worth while 
to plot the point 
corresponding to 


es. 


3 1 Fe) EO RST PORES pee ea Opec 


(iv) When x is any 
number slightly less 
than 1, the sign of 
P'(x)1s 
(<4 +} = (—); 
writing down 
only the signs of the 
various factors. 
When x is any 
number slightly 
greater than 1, the 
sign of P’(x) is 
(pt = 
P'(x) changes sign at 
x = 1, and therefore 
P has a local 
minimum at 1. 
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(v) When x is slightly 


less than — 2, the 
sign of P’(x) is 
(—}(+ }i—} = i+} 
When x is slightly 
greater than —2, the 
sign of P’(x) is 
(—)(+)(—) = (+). 
Hence P has neither 
a local maximum nor 
a local minimum at 
— 2. 


(vi) When x is slightly 
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less than —, the 
sign of P’(x) is 
(—)(+)(—) = (+). 
When x is slightly 
greater than —4, the 
sign of P’(x) is 
(—)(+)(+) =(-). 
Hence P has a local 
maximum at —3. 
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We can now be fairly P(x) 
confident that we can 

join the dotted line to 

complete the sketch. 


Further Exercises 


You will find exercises on this topic in every elementary book on calculus. 
In particular, there are a number of suitable exercises in Apostol, Cal- 
culus Vol. 1 (see Bibliography), if you feel that you need more practice. 


Exercise 1 Exercise 1 
(3 minutes) 


Find the stationary points of the following functions and classify each 
of them as a local maximum, a local minimum or neither. 


(i) g:x+—> 3x* — 4x? (x € R) 
Gi) 5 2+ an x (x €[—7,7)). * 
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Solution 1 


(i) e'(x) = 12(x3 — x?) = 12x?(x — 1), 


g(x) = 01f x = 0 or x =" 1, 
so these are the stationary points. | 


g(x) = 12(3x? — 2x) = 12x(3x — 2), 


SO 
2"(0) = 0 and e”(1) = 12. 

Immediately we can see that there is a local minimum at 1, but Method 

Two breaks down at 0 and we can obtain no information from it. However, 


our Method One will still work. We can see that g’(x) can never be positive 
if x < 1, and so it is certainly negative for x near to and on either side of 0. 


It follows that g can have neither a local maximum nor a local minimum | 


at 0. 
(11) S'(x) = cos x, 


and 
S(x) = Oifx = +5, 


so these are the stationary points, 


S"(x) = —sin x, 
SO 
Tt Tt 
— Beas =] §"|— —- —| 
3)» 


: 7 a 
Hence we can deduce that S has a local maximum at 5 and a local mini- 


: 7 
mam A — 
2 
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Solution 1 


15.2 GEOMETRIC IDEAS 
15.2.0 Introduction 


We shall continue to investigate the problem of optimization, but from 
now on we shall concentrate on functions of two (real) variables; that 
is, functions of the form: 


Fj 2 Ie R= 


where ze R. Before doing this we need to know a little three-dimensional 
co-ordinate geometry, because it is often helpful to represent such func- 
tions by surfaces. We really don’t have enough time to do justice to 
geometry here; in fact a dedicated geometer would probably say that we 
were hardly doing geometry anyway. Our purpose in this section is to 
enable you to visualize the functions, and to describe the corresponding 
techniques in a pictorial and intuitive fashion. In section 15.3 we shall 
apply these geometric notions to our problem of optimization. 


We begin by comparing a problem of optimizing a function of one variable 
with that of optimizing a function of two variables, just to give you a 
concrete example to cling to if you feel that things are getting too abstract. 
Then we shall see how various functions can be represented by surfaces. 
For our purposes, the most important of these surfaces is the plane, and 
the fact that we are able to use planes to investigate more complicated 
surfaces is crucial to the discussion which follows. 


Example 1 


Suppose that an engineer is designing an aqueduct, and he has to use 
rectangular sheets of metal 4 ft wide by 10 ft long. 


10 ft 
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Example 1 


He intends to bend two edges at right-angles, as in the diagram, and his 
problem is to choose the value of x (in feet) which will allow the maximum 
amount of water to travel along the aqueduct. 


We shall try to find the value of the depth x which will give the greatest 
cross-sectional area, A (measured in square feet), of the channel. (It isn’t 
obvious that this is the best thing to do; after all, the shape of the channel 
might alter the speed of the water flowing along it.) This merely requires a 
simple application of the techniques which we have developed in section 
15.1. 

Our problem in mathematical terms is to find the overall maximum value 
of the function: 


j= A (x € [0, 2)). 
We have 
f(x) = x(4 — 2x), 


and therefore 


I(x) = Al — x), 
so that 1 is a stationary point. 
f(x) an <a 


and therefore f has a local maximum at x = 1. It follows that the greatest 
possible value of A is 


A = f(l)=2. 


The above solution is all very well, but wouldn’t the engineer have been 
wiser to bend the sides at an angle other than a right angle? 


In this case it shouldn’t take you long to show that the cross-sectional 
area is (4 — 2x + y),/x? — y’. This time our mathematical problem is to 
find the greatest value of F(x, y), where the function F is defined by: 


F :(x, y)*—> (4 — 2x + y),/x? — y? (x, y)e [0, 2] x [0,2] and y < x). 


The engineer now has two variables to juggle with, x and y. How can he 
solve a problem of this kind? He could use a common-sense approach, 
as we suggested initially for functions of one variable. For instance, he 
could take 10 values of y in all possible combinations with 10 values of x, 
and calculate the value of F corresponding to each pair; but this would 
mean 100 calculations, and if he took 50 values of each variable he would 
require 50* = 2 500 calculations. 


In a practical situation we might have a dozen variables, instead of just 
two, in which case taking only ten values for each would give rise to 
10/17 = 1000000000000 calculations, which would be a formidable 
task for even the fastest modern computer. Someone once said that a 
function of one variable can be tabulated on a page; a function of two 
variables needs a book; a function of three variables requires a shelf of 


34 


FM 15.2.0 


books; a function of four variables would take a whole library; and 
nobody in his right mind would try to tabulate a function of four or more 
variables. 


Clearly we need something equivalent to the technique “a we have 
just developed for functions of one variable. 


We shall restrict the discussion in this unit to functions of two variables, 

although the basic results have equivalent forms for any number of 

variables. (A function of n real variables is a function which maps an 

element of the form (x, y,..., w) to a real number; that is, its domain is a 
ee eee 


nterms $; 


subset of the Cartesian product, Rx Rx.--- x R (which is usually 
nee | 
nterms , 
denoted by R’”), er its codomain i is R.) We shall discuss only functions of 
two variables, because such functions can be visualized in geometric 
terms. 


15.2.1 Graphs of Functions of Two Variables 


We can represent many functions of one (real) variable by (pictorial) 
graphs, which enable us to use our intuition when examining the beha- 
viour of the functions. In particular, when thinking of maxima and minima 
of such functions, we find the graphical approach very helpful. You 
should notice, however, that we try to discard the purely pictorial argu- 
ments in favour of symbolic reasoning, as soon as we feel that we are on 
the right track. 


We shall base a number of arguments on pictures because we think that 
they are easier to understand this way: in later mathematics courses we 
shall need to examine more closely the difficulties which can arise. 


Our first thought is to find a diagram which represents a function of two 
(real) variables, rather as a graph represents a function of one (real) 
variable. In this sort of diagram, we shall find that a function can often 
be represented by a surface. All our functions will be assumed to be 
“well-behaved”; in other words, the surfaces representing them have 
no spikes, gaps, or similar oddities. 


Cartesian Co-ordinates 


In the Cartesian plane each ordered pair of real numbers (a, b) corresponds 
to a unique point. 
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Likewise in the Cartesian space of three dimensions, each ordered triple 
(x, y, z) of real numbers corresponds to a unique point of the space. 


For example, we reach the point (3, 1, 5) if we start at the origin and proceed 


3 units along the x-axis, 
1 unit parallel to the y-axis, 


and 


5 units parallel to the z-axis. 
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Representation of Functions 


For functions of one variable, we know that a function gives rise toa graph Discussion 
(in the sense ofa list), and this can be illustrated as a graph (in the sense of a =e 
picture). We now try to do something similar for a function of two vari- 

ables. Let us look at a particular example of a function of two variables, 

and see how it gives rise first to a list and then to a picture. 


Example 1 Example 1 


Consider the function 
SS eo Oe eS ((x, y)e R x R). 


The ordered pair (3, 4) is mapped to ./37 + 47 = 5, and this corresponds 
to the ordered pair ((3, 4), 5). (Notice that the first element of this pair is 
also a pair.) Similarly, the pair (5, 12) maps to 13, and this corresponds to 
the pair ((5, 12), 13). If we put F(x, y) = z, then (x, y) maps to z, which 
gives rise to the pair ((x, y), z). In this way we can build up a table: 


(x, y) Z 
(3, 4) 5 


(5, 12) 13 


With the pair ((3, 4), 5) we can associate the point with co-ordinates (3, 4, 5); 
with the pair ((5, 12), 13) we can associate the point with co-ordinates 
(5, 12, 13); and so on. In this way, the function defines a set of ordered 
triples. Alternatively, we can think of the equation z = F(x, y) as defining 
a restriction on the variables x, y and z. This restriction corresponds to a 
subset of R x R x R (the set of all ordered triples of real numbers), 
namely the subset {(x, y, z):z = F(x, y)}. 

The surface corresponding to this function F is particularly easy to 
visualize, for in this case 


z = F(x, y) = x* + y? 


Zz (x, y,Z) 


a 


a7 


If we fix the value of z (corresponding to the vertical height in the diagram) 
and look at all the points at this height whose co-ordinates satisfy 


then we find that they are all the same distance, ./x? + y?, from the 
z-axis; that is, they lie on a circle. 


We can describe the surface in words by saying that we move from the 
origin in any horizontal direction, then vertically through the same 
distance to reach the surface. This surface is a cone with its vertex at the 


origin. 
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In general, let F be any function of two (real) variables : 


we, eee 3 ((x, y)e R x R). 


Then to each ordered pair (a,b) we can associate the point Q with co- 
ordinates (a, b, F(a, b)). 


In general, then, a function F of two variables defines a subset of 
R x R x R which is often a surface. It is this surface which is the general- 
ization of the pictorial graph of a function of one real variable. 


Exercise 1 Exercise 1 
(2 minutes) 


Indicate on a diagram the sets of points with co-ordinates (x, y, z) satis- 


fying: 
(i) x = 0, 
(ii) y = 0, 
(x = -7 = 6: = 
Exercise 2 Exercise 2 


(2 minutes) 
Mark on a diagram the set of points with co-ordinates (x, y, z) in 
R x R x R which corresponds to the condition: 


2x —y=0. = 
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Solution 1 Solution 1 


Solution 2 Solution 2 


Considering the equation as a restriction defining a subset of R x R x R, 
we see that z does not appear in the equation, so there is no restriction on 
z. But x and y are restricted. If (x, y, z) is to belong to the subset, then x and 
y must satisfy the equation 2x — y = 0. The set of triples (x, y, 0) which 
satisfy this equation from a line in the xy-plane. 


Corresponding to any point on this line, we can get other elements of the 
required subset of R x R x R by choosing any value of z. All in all, 
we get a plane perpendicular to the xy-plane, which intersects the xy- 
plane in the line with equation 2x — y = 0. 
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The plane y = 2x 


15.2.2 The Intersection of Surfaces 15.2.2 


Before we find the equation which defines a general plane, we would like Discussion 
to give you some reasons for our interest in the subject. At first sight the EE 
following examples have nothing to do with planes, but a closer examina- 

tion will reveal the connection. — 


Contour Lines 


It is difficult to give an impression of a three-dimensional object on a two- 
dimensional piece of paper. One way of overcoming this difficulty, which 
is shown in the following examples, may be quite familiar to you, and we 
can develop it into a useful tool for examining functions of two variables. 


Example 1 Barometric Pressure Example 1 


ATLANTIC WEATHER MAP 


SITUATION AT NOON JAN.22 


The red curves on the weather map join the points of equal barometric 
pressure. (The pressure varies with height, but the values shown refer to 
the pressure at sea-level.) The function illustrated in this case is 


P :(point on the map)+—— (barometric pressure at the corre- 
sponding point on the earth’s surface). 
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Example 2. Ordnance Survey Maps Example 2 


Zs 
LOCKERBROOK 
HEIGHTS 
= 


7S 


The cartographer has only a flat piece of paper, but he does his best to 
give an impression of the shape of the land surface by showing us the 
contour lines; in other words, he joins the points of equal height above 


sea level. 


The function which is illustrated in this case is 


Example 3 


h:(point on the map)+—— (height of the corresponding point 


above sea level) 


Example 3 


Suppose that we take a circular cylinder of unit radius which has its axis 
vertical (along the z-axis) and intersect it with a horizontal plane 5 units 
above the xy-plane, as shown in the following diagram. 
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The two surfaces (the cylinder and the plane) meet in a curve, which is in 
fact a horizontal circle 5 units above the xy-plane. 


Any point on the cylinder is one unit from the z-axis and therefore 


Jey = |. The equation of the cylinder is therefore ./x? + y? = 1, 
by which we mean that the set of all points with co-ordinates (x, y, z) 
inR x R x R = R®, satisfying this equation, lie on the cylinder. 


The equation of the plane is z = 5, and the two equations taken together : 
af a = 
z= 5 


determine the set of points lying on the red circle. 
Another way of writing this is as follows: denote a point P with co- 
ordinates (x, y, z) by P(x, y, z); then the cylinder is the set 


A = {P(x, y, 2):(x, y,z)E R*,,/x* + y = 1}; 
the plane is the set 
B= Pee Be er c= 5: 
and the circle is A 7 B. 
a ; 
You may say that our two equations imply that x* + y? = 5° so isn’t this 


the equation of the circle? If so, compare the following two sets in 
Kx Kx KE: 
(i) A © B, the set of points with co-ordinates (x, y, z) which satisfy 


/x+y=1 
and 


a= 3- 
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(ii) C, the set of points with co-ordinates (x, y, z) which satisfy 


2 = 
ey => 
Notice that every point in AB has z = 5, so that it must have co- 
ordinates of the form (x, y, 5), whereas the point with co-ordinates (2, 2, 40), 
for instance, lies in C but not in A B. Every point in Ao B satisfies 
the conditions for C, but the converse is not true, so that A 4 Bis a proper 
subset of C. This is to be expected because A B implies a stronger 
restriction in R x R x R than does C. To belong to C, (x, y, z) must 


be such that x? + y? is the same as z To belong to A - B, (x, y, z) must 


‘ Z ; 
be such that, not only is x* + y* the same as 5° but also both expressions 


have the value 1. In fact the points in C lie on what is called a paraboloid 
of revolution, and this paraboloid contains the red circle. 


A Generalization Discussion 


Suppose that we are given an arbitrary function of two variables, F, with 
domain the xy-plane, and we intersect the surface 


z= (x, y) 
with the horizontal plane 


Z=C 
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v 
y 
x 
Various contour lines shown in red for various values of c. 
The resulting curve is called the contour line corresponding to the height Definition 1 


* 


c. Taking various values of c will give a set of contour lines, which, when 
viewed from above (looking down the z-axis), could look like this: 


This is simply a general version of the pressure and altitude diagrams which 
we used to introduce these ideas. 


The previous examples have shown how planes parallel to the xy-plane 
(horizontal planes) can be used to describe surfaces, but we intend to use 
planes parallel to the z-axis (vertical planes) too. 


Consider the intersection of our arbitrary surface defined by 
g = Fy) 
with the plane 


y=2. 
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The effect is rather like slicing a Dutch cheese; for each slice the red rind 
of the cheese forms a different curve, and for each value of b in the above 
diagram we get a new red curve. The advantage of this idea is that it 
reduces a surface, which is difficult to draw, to a set of curves, each lying 
in a plane, which we can draw on a piece of paper. Mathematically speak- 
ing, we have reduced a function of two variables to a whole set of functions 


of one variable each corresponding to a particular value of b, and a particu- 
lar red curve. 
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15.2.3. The General Equation of a Plane 15.2.3 


Up to this point we have been finding the surfaces corresponding to given Discussion 
equations and functions, but now we want to put the problem in reverse. _ 

Can we find equations of given surfaces? What is the equation of a plane? 

If you have already seen the television programme, then you will know 

that this is an essential step on our way to solving optimization problems 

for functions of two variables. 


Suppose that the plane passes through the point (a,b,c) and that it is 
inclined at an angle A in the x direction and an angle B in the y direction. 
You can see what this means from the cut-out diagram provided,* which, 
when set up with a sheet of paper placed on top to represent the plane, 
looks like this: : 


et oe ae 


The value of z corresponding to an arbitrary choice of (x, y) is simply 
the result of adding the three terms on the right of the above diagram, 


z=c+(x—a)tanA +(y — b)tanB, Equation (1) 


which is the required equation of the plane if neither A nor B is > 


If either A or B is > then the plane is perpendicular to the xy-plane, 


and looking down the z-axis we see only the line PQ. 


* The cut-out diagram is included inside the back cover of this text. 
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_ In the xy-plane we could represent the line PQ by the equation 
ax + Py+6=0. | 


In R x R x R this equation represents the plane. 
Equation (1) and Equation (2) are particular cases of the equation 


ax+ By +yz+6=0 


where «, B, y and 6 are real numbers independent of x, y, z (see also the 
following exercise). This is the general equation of a plane. Notice par- 
ticularly that the plane is horizontal (that is, parallel to the xy-plane) if 
a = B = 0(in other words when the angles A and B are both zero). 


Exercise 1 


(i) What values should we take for a, f, y and 6 in order to make Equation 
(3) identical to Equation (1)? 

(ii) At what points does the plane ax + By + yz + 6 = 0 meet each of 
the three co-ordinate axes? 

(iii) The equation A(x — a) + p(y — b) = O represents a plane perpendic- 
ular to the xy-plane which passes through the points (a, b, z) for any 
value of z. What effect does it have on the plane if we vary the values 
of A and pn? 


(A and ware the Greek letters called “lambda” and ‘‘mu” respectively.) 
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Equation (2) 


3 Equation (3) 


Exercise 1 
(3 minutes) 


15.3 PARTIAL DERIVATIVES AND OPTIMIZATION 
OF FUNCTIONS OF TWO VARIABLES 


15.3.0 Introduction 


We have derived the general equation of a plane, but if you have seen the 
television programme, you will know that we really need the equation 
of the tangent plane at a point on a given surface. We can then imagine 
this plane moving over the surface, and we hope that this notion will give 
us a technique for finding the maximum (or minimum) value of the 
corresponding function, just as a similar idea helped for functions of 
one variable. For the moment we need something like the derivative of a 
function of one variable, which was useful when discussing rate of change. 
The corresponding concept which we are going to examine is that of a 
partial derivative. 


First let us give an intuitive idea of the concept of partial derivative. 
Imagine yourself standing at a crossroads on a hillside, the roads running 
East-West and North-South. Roughly speaking, the slopes of the East— 
West road and the North-South road are the partial derivatives of the 
function, represented by the hillside, at the point where the roads cross. 
If the crossroads happened to be at the top of a hill then each of the slopes 
would be zero. It is this intuitive idea that we want to make precise, and 
the following example will lead us in the right direction. 


The geometric discussion in the example is intended to help you to 
understand the definitions which follow. Some people find three-dimen- 
sional figures hard to visualize, and if you don’t like the geometry, you 
may be better off going directly to the definitions. You certainly should 
not spend a great deal of time trying to understand it if you find it difficult. 


Example 1 
Consider the surface representing the function 


F:(x,y)—> f/1-(@? +?) (uy)eR xX Rx’? + y’ <1). 


The domain of F is represented in the xy-plane by the points on and within 
the circle with radius 1, centred at the origin. 


If we let z = ./1 —(x? + y?), then it follows that x? + y? + 27 =1. 
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15.3.0 


Introduction 


x * 


Example | 


Equation (1) 


(continued on page 50) 


Solution 15.2.3.] 


(i) a =tan A 
6 =tan B 
7 =) 


6 =c—atanA — btanB. 
(11) If « # 0, the plane meets the x-axis in the point with co-ordinates 


[—2.0.0] 
a 


If « = 0, the equation of the plane is By + yz + 6 = 0: 

(a) If 6 # 0, then y and z cannot be zero simultaneously, so the plane 
does not meet the x-axis. 

(b) If 5 = 0, the equation of the plane is By + yz = 0. All points with 
co-ordinates (x, 0, 0) lie on this plane, and so the plane contains 
the whole x-axis. 

A similar argument can be used for the other co-ordinate axes. 

(iii) Varying / and w alters the line in which the plane cuts the xy-plane. 

It will not change the fact that the plane is perpendicular to the Xy- 


plane. = 


(continued from page 49) 


_ The distance of any point, P(x, y,z), in R x R x R from the origin is 
= x? + y? + z?. This can be seen in the diagram, first by using Pythag- 
oras’s Theorem in the triangle ONM, and then in the triangle OMP. 
Since points on the surface satisfy the equation x? + y? + z? = 1, it 
follows that any point P lying on the surface must be at unit distance 
from the origin, and since z is always positive, Equation (1) represents a 
hemisphere. 

In terms of our intuitive discussion, this hemisphere is the hillside. We 
are now going to choose a point Q on this hillside and assume that Q 
is our crossroads, with the roads through Q in the planes parallel to the x 
and y axes. Any point Q will do, and to illustrate the idea we choose QO to be 


. ; 1 
the point with co-ordinates | 5° 0, — 
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1 3 
The point Q with co-ordinates > 0, ve lies on the surface and on the 


plane y = 0. If we were to cut the hemisphere with the plane y = 0 
through Q, and then look along the y-axis, we would see the semi-circle 
shown in red in the following diagram: 


The plane y = 0 


The plane y = 0 


You can imagine the red curve to be the road through Q in the x-direction ; 
later we shall find it very useful to be able to calculate the slope of such 
curves at an arbitrary point Q. 


Next we find the slope of the red semi-circle at the point Q with co- 


1 3 = 
ordinates > 0, “7 |" The semi-circle is determined by the equations : 


gagi= +7) 
y=0 
so that on the curve we have: 
z= Ji =x" (xe[—1, +1)). 


But a relationship like this defines a function, f, say: 


fox /1i - x (xe[—1, +1)). 
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Since we are interested in the slope at Q, we calculate the derivative, f(x). 
We get 


Fi) = 


l1—x 
which takes the value hen x 
——=w =, 
J3 2 
Suppose now that we want the a = the other road at Q. First intersect 
the hemisphere with the plane x = 3 


F The plane 


The red curve in the diagram is determined by the equations: 


r= JI) 


xX => 


dOo|-— 


so that on the curve we have: 
PF (of) 
This relationship defines a function, f, say: 
fe sees = [ye 3 2}. 
Again we calculate the derivative, f(y). We get 


f= 


ae 
which takes the value 0 when y = 0 (the slope of the curve at Q seemed 
likely to be zero from the diagram, so this should be no surprise). 
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15.3.1 Definition of Partial Derivatives 


Having worked through the previous example, you may well feel that there 
must be a quicker way, and indeed there is. It is, however, the intuitive 
ideas of that example which point the way. 


Suppose that we examine the function F from a slightly different point of 
view. We know that 


F :(x, y)—> ./1 — (x? + 2) (x y)e (x y):x? + y* < 1}). 


If we keep y constant, y = b say, then we obtain a new function (of one 
variable): 


fix /1 —(@? +b) (xe[—-/1 — b?, + \/1 — Bb’), 
and 
—x 
52S ares 


Similarly, if we keep x constant, x = a say, then we obtain a new function 
(of one variable): 


fairy —> 1 -(@ + yy?) (yve[-J/1 - a’, + /1 — a’), 


and 


1 3 
The slopes of the roads through the point Q, e 0, 4 running parallel 


to the x and y axes are given by the derivatives f \(5) and f5(0) respectively 
(taking a = 5 and b = 0). 
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15.3.1 


Definitions 


*-*-s 


The expression f(x) gives the slope of the surface (defined by the func- 
tion F) in the direction of the x-axis at the point (x, b): that is, F(x) is the 
rate of change of F with respect to x, when y has the constant value b. 


Similarly, f5(y) gives the slope of the surface in the direction of the y-axis 
at the point (a, y); that is, f(y) is the rate of change of F with respect to y, 
when x has the constant value a. 


1 3 
We chose to consider the point Q with co-ordinates | > 0, I SO we took 


a = 4 and b = 0. We would like to know the corresponding slopes (rates 
of change) at any point on the surface defined by F; that is, we now wish 
to vary a and b. This means that we need to express the slopes in terms of 
functions of two variables. So we define two new functions, F, and F;, 
by the equations: 


FS Ope eas ((x, y) © {(x, y):x? + y? < 1}) 
and 
(x, y) = = ((x, y) € {(x, y):x? + y? < 1). 


1 — (x? + y”) 


We are thus led to the following definition of partial derivatives of a 
function, F, of two variables x and y. 


The partial derivative of F with respect to the first variable, x, at (x, y) 
1S 


F(x + h, y) — F(x, y) 


F’,(x, y) = lim 
h~0 h 


The partial derivative of F with respect to the second variable, y, at (x, y) 
is 


- ° F(x, +k} = #6, ) 
F(x, y) = lim = 


In order to find the two partial derivatives, we simply keep each of the 
variables fixed in turn and differentiate with respect to the other. 


Example 1 
If | 
G5 9) ay FX ((x, y)e R x R), 


then, differentiating with respect to x at (x, y), we regard y as constant and 
get 


G'(x, y) = 2y + 2x, 


and differentiating with respect to y at (x, y), we regard x as constant and 
get 


a(x, y) = 2x. 
(Working directly from the definitions, we have: 


2Wx + h) + (x + h)*? — (2xy + x?) 
h 


2yh + 2xh + “| 


G(x, y) = lim | 


x: fam 
h~0O 


h 


= 2y + 2x, 
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Definition I 


ae ee 


Definition 2 


kek 


Example I 


and 
+ —— 2 2 
ie the x(y+k)+x (2xy + x”) 
k~0 k 
ORE posse 
at = * 
= 2x = 
Exercise 1 


Find the partial derivatives at (x, y) of the functions defined by the follow- 
ing equations; each function has domain R x R. 
(i) F(x,y)=x?+y’ 
(11) G(x, y) = x exp (xy) 
(ii) H(x, y) = x sin(x + y) 
(vy) Pe 7) = x* 4 -y* — Axt*y”. 2 


Alternative Notation 


There are various notations for the partial derivatives ; the most common 


_ OF : ; 
is = for what we write as F(x, y). This notation arose presumably because 
x 


——- —— 
of the commonly used notation £ for the derivative of a function, f, 
x 


OF 
of one variable. If you use the _ notation, then you must be extremely 
x 


5 oF 2 = 
careful later. For example, it is not generally true that ra ; = is simply 
* 


F : 
op as the notation would suggest. 


The alternative notation F,, or F, is also very common. 
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Exercise 1 
(2 minutes for each part) 
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Solution 1 Solution 1 
(i) F4(x, y) = 2x 
F(x, y) = 2y 


(1) G(x, y) = exp (xy) + xy exp (xy) 
2(X, y) = x* exp (xy) | 

(111) H4(x, y) = sin(x + y) + xcos(x + y) 
(x, y) = x cos(x + y) 

(iv) P(x, y) = 4x3 — 8xy° 


P5(x, y} = 4y" = 12x" 9: = 
15.3.2 The Tangent Plane 3 15.3.2 
You probably have an intuitive notion of what we mean by the tangent Discussion 


plane at a particular point on a surface. It is, after all, the plane which sits 
comfortably on the surface at the point in question. Once again, we 
assume that our surfaces are smooth with no sharp projections. It would, 
for example, be difficult to decide where the tangent plane should be on 
the apex of a church steeple. 


On the other hand, it is quite easy to imagine a tangent plane at a point 
on a smooth sphere. 


56 


We shall now define the tangent plane at any point on a smooth surface. 
Suppose that we are given a surface defined by: 


F :(x, y)+—> F(x,y) — (x, ye R x R), 


and we wish to define the tangent plane at the point Q with co-ordinates 
(a, b, F(a, b)). We have seen in Exercise 15.2.3.1 (iii) that the equation 


A(x — a) + ply — b) =0 


defines a plane which passes through Q and is perpendicular to the xy- 
plane. 


cg surtane SF ew) gurtace z=F (x,y) 
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The intersection of this plane with the surface will be a curve (shown in 
red on the diagram). This curve passes through Q and has a tangent line 
(shown by a heavy black line) at Q. If we vary the values of 2 and y, the 
plane will rotate about the vertical line through Q, and each pair of values 
of 2 and u will give us such a tangent line. If all these tangent lines at Q lie 
in a plane, then we call this plane the tangent plane at Q. 


The Equation of the Tangent Plane 


Our assumption that the surfaces we meet are smooth is intended to 
imply that there is a tangent plane to the surface z = F(x, y) at O, but how 
can we find its equation? 


Suppose that we take the particular values uw = 1, A = 0, in the equation 
Ax — a) + u(y — b) = 0. 
Then we simply get the equation 

y=, 
and the slope of the corresponding curve of intersection at Q is F‘,(a, b). 
In other words, F’,(a, b) is the slope of the tangent to this curve at Q. 


The curve of intersection is shown in red on the following diagram, and 
the tangent to this curve at Q is shown by a heavy black line. 
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Definition 1 


xk 


Main Text 


x* * 
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The plane y=b 


“~The slope of this 
line is F;, (a,b) 


ere or ne ame ap Ht EY A NN ene art be be See 


XN 


The slope of the tangent line in the above diagram is simply F(a, b), and 
similarly the slope of the corresponding line in the y-direction is F(a, b). 


The slope of this 
line is F, (a,b) 


"~The plane x=a 


ae a ae ae ae et me 


We can now use our cut-out model of a plane again, this time to find the 
equation of the tangent plane at Q. 


- The slope of this 
line is F; (a,b) 


The slope of this 
line is F, (a,b) 
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In the formula on page 47, we simply put tan A = F(a, b), tan B = F(a, b) 
and c = F(a, b). We can see from the last diagram that any point on the 
tangent plane has co-ordinates (x, y, z) which satisfy the equation 


z = F(a,b) + Fy (a, b)(x — a) + F%(a, b)(y — b) 


and this is the equation of the tangent plane to the surface at (a, b, F(a, b)) 


Exercise 1 


For each of the following functions, find the equation of the tangent plane 
at the point on the surface corresponding to the pair (a, b) (each function 
has domain R x R). 
Brae x + 
(11) G :(x, y)+—— x exp (xy) 
(ii) H :(x, y)-——> x sin (x + y) 
(iv) P:(x, y)——> x* + y* — 4x’y”. 
(You will be able to use the results of Exercise 15.3.1.1.) s 


Exercise 2 


Find the equation of the tangent plane to the surface defined by 
F :(x, y)->- (4 — 2x + y)./x? — y’ 
((x, y)€ {(x, y):0 < y < x < 2}) 


at the point (a, b, F(a, b)). (This is closely related to the aqueduct problem 
which we discussed in section 15.2.0.) & 
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Equation of 
Tangent Plane 


kk 


Exercise 1 
(2 minutes for each part) 


Exercise 2 
(3 minutes) 


PM -I3.3.2-15.35 


Solution 1 ~ Solution 1 


(i) z = a? + b* + 2a(x — a) + 2b(y — BD) , 
(ii) z = aexp (ab) + (ab exp (ab) + exp (ab))(x — a) + a? exp (ab)(y — b) 
(ill) z = asin(a + b) + (sin(a + b) + acos(a + b))(x — a) | 
+ acos(a + b)(y — b) 
(iv) z = a* + b* — 4a’b? + (4a? — 8ab*)(x — a) 
+ (4b? — 12a*b?)(y — b). = 


Solution 2 , | Solution 2 


First, we must find the partial derivatives; they are: 


xX 
10% y) = —2/x? = y? + (4 = 2x + YH 
ee: 


a. 


: 
Fox, y) = /x? — y? — (4 - 2x + J 
af XSF 


Thus the equation of the tangent plane is 


z = (4 — 2a + b)./a* — b? 


Ss 
+ = 2/8 + Se (x — a) 
| 4—2 b 
tiiga = a 


Je =F 


Unfortunately, this cannot be simplified ; we included the exercise because 


of its tie-up with the aqueduct problem. = 
15.3.3. Optimizing Functions of Two Variables 15.3.3 
There is no doubt that finding the overall maximum (or minimum) of a Discussion 


x* 


function of two variables is, in general, harder than finding the overall 
maximum (or minimum) of a function of one variable. Some might say, 
“more than twice as hard’’. For a differentiable function f of one variable 
with domain [a,b] we merely locate the stationary points, investigate 
their nature, and then examine the values of f(a) and f(b). The domain of 
a (real) function of two variables is a subset of R x R, and instead of just 
two end-points we are now likely to have a curve as the boundary of our 
domain. 
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For example, we have already discussed the function 
Pix yy— /1- 6? +) i yet + ST) 


The domain has the circle x* + y? = 1 in the xy-plane as its boundary. 


Suppose that we wish to find the overall maximum (or minimum) value 
of the images ofa function F with domain, A, a subset of R x R. The points 
where the tangent plane is parallel to the xy-plane, on the surface defined 
by z = F(x, y), are clearly going to be of interest. This leads us to our 
next definition. 


If F(a, b) = 0 and F(a, b) = 0, then (a, b) is called a stationary point of F. 


Notice that since the tangent plane at a stationary point 1s parallel to the 
xy-plane, its equation is simply z = F(a, b). 


15.3.4 Local Maxima and Minima 


You may find the precise definitions of local maximum and local minimum 
a little hard to digest, so we give intuitive definitions first. 

If (a, b) is a point in the domain of F, and if F(x, y) < F(a, 5) for all (x, y) 
in the domain of F sufficiently close to (a, b), then we say that F has a 
local maximum at (a, b). 

If (a, b) is a point in the domain of F, and if F(x, y) > F(a, b) for all (x, y) 
in the domain of F sufficiently close to (a, b), then we say that F has a 
local minimum at (a, 5). 


Speaking very roughly, if it rains on the surface, the puddles collect 
around the local minima, and the water runs away to the overall minimum 
as the puddles overflow. 


The difficulty with the above definitions is that they depend on the meaning 
of “‘sufficiently close’, and it is this phrase which needs to be precisely 
defined. 


If we use our approach to functions of one variable as a guide, then we 
need a “‘small” set in R x R where before we had a “‘small’’ interval, 
fc — e,c + €] in R; the most suitable set in R x R is a circular disc. 

We let S(a, b, «) denote the set {(x, y):(x — a)? + (y — b)? < &*}, which is 
a disc with centre at the point (a,b) and radius ¢. 


Let F be a function with domain A CR x R. Then, following our 
definitions for functions of one variable, we make the following formal 
definitions. 
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Definition 1 


15.3.4 


Definitions 
x“ xr 


Intuitive Definition 1 


xk*k 


Intuitive Definition 2 


xk 


Discussion 
x 


If there is a positive number ¢« such that F(x, y) < F(a,b) for all 
(x, y)e AO S(a, b, €), then we say that F has a local maximum at (a, b). 


If there is a positive number e such that F(x, y) > F(a,b) for all 
(x, y)€ A S(a, b, €), then we say that F has a local minimum at (a, b). 


Let us have a look at the set A 7 S(a, b, e). Take the set in the following 
diagram to be the set A. 


If we represent S by a 10p coin, then, placing the 10p down on the set, the 
part of A which is covered by the coin is the set A A S(a, b, e). The point 
(a, b) is of course the centre point of the coin, and ¢ is its radius. 


Example 1 
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(ii) 


Local maximum 


2 a 


An S(a,b,e) 


| 


z 
Overall maximum 
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Definition |! 


ce MB 


Definition 2 


x*r* 


Example I 
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Just as for functions of one variable, we have two major problems. A Discussion 
stationary point need not be a local maximum nor a local minimum. 
Stationary points of this kind are called saddle points (for obvious reasons). Definition 3 


Horizontal tangent plane at a saddle point 


In other words, saying that the tangent plane is horizontal guarantees 
neither a local maximum nor a local minimum. 

The second big problem is that a local maximum, or indeed an overall 
maximum, can occur on the boundary of the domain, and similarly for 
local and overall minima. If we restrict ourselves to a search for stationary 
points, then we may miss points of this kind. 


We shall attempt to overcome the first problem, but we shall not have 
time to find a way to overcome the second problem in this course. In 
general, the second problem is not so serious when applying the methods 
to physical situations, since the information from the situation itself 
often helps us to decide the nature of the point we are dealing with. 


15.3.5 The Aqueduct Problem 15.3.5 


We are at last in a position to attempt to solve the engineer’s problem of Example 
designing the aqueduct which we mentioned on page 33. We left him : 
trying to decide on the choice of x and y which would give the overall 

maximum value of the function: 


F :(x, y)-—> (4 — 2x + y),/x? — y*, (x, y)€ [0,2] x [0,2] and y < x). 


It is important to notice that we do not need to draw a picture of the surface 
defined by the equation 


z= (4 — 2x + y),/x? — y’, 


in fact, we want to avoid this if possible. 
We found in Exercise 15.3.2.2 that 


(a,b) = —2,/a? — b? + Faa)4 - 8 + 


and 
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The tangent plane is horizontal if F(a, b) = F(a, b) = 0, so we need to 
solve the pair of equations: 


=D =| ial 2B =) 
qe 


which simplifies to 


(4— 2a + b) =a? — b’ 


No} s 


b4 — 2a + b) = a? — b?) 


a : 
These equations are satisfied simultaneously if b = > and then, sub- 


stituting for b in the first equation, we obtain a = 4, and therefore b = z. 


(The other* solution, a = b = 0, clearly corresponds to a minimum value 
of F(x, y), namely 0.) The corresponding value of the cross-sectional area is 


4 —— 
aE ~ 2.3, which is obtained by replacing x by $ and y by 3 in the expression 


for F(x, y). (Remember that when the sides of the aqueduct are vertical 
the greatest value of the cross-sectional area is 2.) 


The engineer would surely find this result very impressive. Just by doing Discussion 
a fairly simple calculation, he would be able to increase the flow of water 

by about 15% at no extra cost in material. He has a greater value for the 

area than before, so there is some cause for satisfaction. It would, of course, 

be even better if he could assure himself that he had found the greatest 

possible value. 


In this case we are fairly certain from the physical origins of the problem 
that the point where the tangent plane is horizontal will give a local 
maximum value of F. However, we cannot always be this certain, and it 
would help to have a technique which would distinguish between local 
maxima, local minima, and other points where the tangent plane is 
horizontal. This is our next task, and we shall begin by considering a very 
easy problem for which we know the answer in advance. 


15.3.6. A Useful Technique 15.3.6 


Example 1 | Example 1 
Find the overall minimum value of the function 


F :(x, y)-——> x2. + 2y*_ (x, ye R x RI. 2 


Solution of Example 1 


Since both expressions on the right are positive or zero, the answer is 
obviously 0, and this occurs when x = y = 0. 


Now let us use the above example to test a technique for classifying Discussion 
stationary points, which we can use when the answer isn’t obvious. =e 


* The solution a = b = 4 does not correspond to a point in the domain of F. 
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The surface z = x* + 2y? looks like this: 


Notice that 
F(a, b) = 2a and _ F(a, b) = 4b, 


and therefore the point corresponding to a = b = 01s a stationary point, 
confirming what we already know. 


The essential idea of our proposed technique is as follows: intersect the 
surface with the plane whose equation is 7 


y=x tang, 


to give the curve shown in red. 


It is obvious from the diagram that the red curve has a minimum at 0, but 
can we show that this is the case mathematically (thus making the diagram 
redundant)? 
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The red curve is determined geometrically as the intersection of the plane, 
{(x, y):y = xtana}, and the surface, {(x, y):z = x* + 2y7}, or, more 
briefly, it is determined by the equations 


x? + 2y? 


Z 


y= x tana. 


Ifr is the hypotenuse of the right-angled triangle with sides x and y (shown 
on the previous diagram), then the second equation can be replaced by 


x =rcosa and y=rsin«a. 
On the red curve we then have 


z = r*(cos? a + 2sin? a) 


r*(1 + sin «), 


so that on the red curve, z(= F(x, y)) takes its minimum value when r = 0. 


In more difficult cases we might have to adopt the following line of 
reasoning to achieve the required result. 


The equation z = r7(1 + sin’ «) defines a function of one variable: 


gr z FER) 
Differentiating, we obtain 

o'(r) = 2r(1 + sin? «) 
and 

O°) = 20 + aa @. 


In particular, #’(0) = 2(1 + sin? «), so that 6"(0) > 0, and therefore the 
red curve has a local minimum at r = 0. 


The essential point about the equation #"(0) = 2(1 + sin? a) is that it 
shows that #’(0) > 0 for all values of «, so that all possible red curves 
obtained in this way have a local minimum at 0. This seems to show 
pretty conclusively that F has a local minimum at 0 (again confirming 
the result which we already know). = 


Roughly speaking, the technique can be summarized as follows. Slice 
through the point of interest on the surface with a vertical cut, revealing a 
curve like the red curve of our example. If all such curves have a local 
minimum at the point, then we would expect the surface to have a local 
minimum there. Our next example will again illustrate this idea. 


This technique for classifying the stationary points will be adequate when 
the surface is “‘smooth’’, and it will certainly be adequate for all the 
problems which you will meet in the Foundation Course. However, it is 
an amazing fact that one can construct a surface for which all the red 
curves have a local minimum at the origin, and yet the surface does not 
have a local minimum at that point. Such a surface cannot be “smooth” 
in our sense, and you may like to try to think of such an example. 


Example 2 (This example has applications in statistics.) 


Given the three points with co-ordinates (1, 1), (2, 3) and (3, 4), find a line 
with equation y = mx +c, such that the sum of the squares of the 
‘vertical’? distances of the points from the line is a minimum. 3 
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Example 2 
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Solution of Example 2 
In the diagram, d,, d, and d, are the vertical distances, and we want to 
minimize 

d?+ 434+ d3 =(m+c-— 1) + (m4 c — 3)? + 3m+4+c — 4)’; 
we can therefore define a function F with domain R x R by putting 

F(m, c) = (m+ c — 1)? + (2m + c — 3)? + 3m +c — 4). Equation (1) 
We then have 

F‘(m, c) = 2m + c — 1) + 42m + c — 3) + 6(3m 4+ c — 4) 

= 28m + 12c — 38 


and 


F4(m, c) = 2m + c — 1) + 22m 4+ c — 3) + 23m +4 c — 4) 


12m + 6c — 16. 


The values of m and c for which F‘,(m, c) = F5(m, c) = 0 are determined 
by the equations 


14m + 6c — 19 =0 
12m + 6c — 16=0 
from which we deduce that 
; and c= —j. 


These values of m and c determine a stationary point of F, but is this 
stationary point the overall minimum? Could we show that the stationary 
point is even a local minimum? Let us try the method of our previous 
example. 


Our variables are m and c instead of x and y, and so we intersect the un- 
known surface by planes perpendicular to the mc-plane and through the 
point (3, — 4). A typical plane has equation 


c+4=(m— 3)tana 
and just as before we can replace this equation by 
m=3+rcosa c=-—i4+4+rsina 


where r is now the distance of the point with co-ordinates (m, c) from the 
point (3, — 4). Substituting these expressions for m and c into Equation (1), 
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we obtain the function ¢ defined by 
P(r) = (r(cos a + sin a) + §)? + (r(2. cos « + sin a) — 3)° 
+ (r(3 cos « + sina) + Z)* rer 3. 
Differentiating twice 
(0) = 2(cos « + sin «)? + (2cos « + sin a)” 
+ (3 cos « + sin a)’) 
and therefore @’(0) > 0 for all values of « (since the squares cannot be 


zero simultaneously); hence the stationary point is a local minimum. 


Since there is only one stationary point and it is a local minimum, it seems 
very likely that we have indeed found the required values of m and c. 
So the equation of the required line is 6y = 9x — 2. 


The difficulty with points on the boundary of the domain of F does not 
occur in this case, because the domain is the whole set R x R, and there 
are no boundary points; but to be safe we ought really to find the images 
of the function when r is very large. For the moment we shall avoid this 


difficulty too. = 

Exercise | | Exercise 1 
: : a (4 minutes) 

Find the equation of the red line which gives the minimum value of 

Pi + p2 + P3. 

Exercise 2 Exercise 2 


(3 minutes) 
Find the point P for which the sum S? + S43 + S% isa minimum. 
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15.4 SUMMARY AND CONCLUSIONS 
15.4.1 Summary 


Our original problem was to find the overall maximum value attained by 
the images ofa function f, This problem is difficult, but we can go some way 
towards a solution if we can locate the local maxima. We know that the 
overall maximum value belongs to the set of local maximum values; it is 
the greatest element of that set. 


The set of local 
maximum values 


@ The overall 
maximum value 


We hope to locate the local maxima by first locating the stationary points 
and then examining their nature, but there are two complications: 


(i) local maxima at the boundary of the domain need not be stationary 
points; 
(11) stationary points can also be local minima or saddle points. 


Local 
Minima 


Local 
Maxima 


Stationary 
Points 


Saddle 
Points 


We can locate the red set fairly easily, using the various methods which 
we have developed. For a function of one variable with an interval [a, b] 
as its domain there are are only two boundary points, and we can use 
the strategy which we developed in section 15.1. For a function of two 
variables the situation is more difficult, for in this case there are, in general, 
an infinite number of boundary points. Our previous strategy will not 
work in this case, because we cannot test an infinite number of values to 
find the largest. But often the physical origins of the problem give us 
useful information which allows us to resolve the problem. | 
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15.4 
15.4.1 


Summary 


Solution 15.3.6.1 


Writing the equation of the line in the form y = mx + c leads to some 
untidy algebra. In this exercise it is convenient to take x = my + cas the 
equation of the line; we then have: 


pit ps t+ p3 =(1—(m4+ 0)? +(2— (3m + 0))? + (3 — (4m + o))? 
= 26m? + 16mc + 3c? — 38m — 12c + 14. 
So F is the function defined by 
F :(m, c)—— 26m? + 16mc + 3c? — 38m — 12c + 14 
((m,c)ER x R)). 
Thus, for a stationary point, we have 
‘(m, c) = 52m + 16c — 38 = 0 
2(m, c) = 16m + 6c — 12 = 0. 


The solution of this pair of equations is 


Nv 


ae =< 
im =a; = 


An argument similar to that given in the text would show that these 
values of m and c do give a local minimum. Thus the equation of the line is 


14x = 9y + 4. = 


Solution 15.3.6.2 
The sum | 
Sj + S3 + 83 =(x — 17 + (y — 2)? + (x — 2)’ 
iy a + Oe 
= 3x? + 3y? — 14x — 12y + 35. 
So we let F be the function 
F :(x, y)-——>- 3x? + 3y? — 14x —12y +35 = (x, ye RR x R). 
We have 
F'(x, y) = 6x — 14 
F(x, y) = 6y — 12. 
Thus for a stationary point, x = 23 and y = 2. 


Once again we can use either the argument given in the text or the geometry 
of the situation to convince ourselves that (24,2) does in fact give the 
point P for which the sum is a minimum. S 
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Solution 15.3.6.1 


Solution 15.3.6.2 


15.4.2 In Conclusion 


Our concern in this unit has been almost entirely with the optimization 

of functions, and our arguments have been based wholly on plausible 

geometric notions. As we have mentioned several times, things can go 

wildly wrong. 

Our reasoning has been almost entirely heuristic, and it is worth quoting 

Polya on this point (see Polya page 113): 
‘‘Heuristic reasoning is reasoning not regarded as final and strict but 
as provisional and plausible only, whose purpose is to discover the 
solution of the present problem. We are often obliged to use heuristic 
reasoning. We shall attain complete certainty when we shall have 
obtained the complete solution, but before obtaining certainty we must 
often be satisfied with a more or less plausible guess. We may need the 
provisional before we attain the final. We need heuristic reasoning 
when we construct a strict proof as we need scaffolding when we erect 
a building... Heuristic reasoning is good in itself. What is bad is to 
mix up heuristic reasoning with rigorous proof. What is worse is to 
sell heuristic reasoning for rigorous proof.”’ 


You may have noticed that we have been careful to say that our results 
are “likely” or “‘plausible’’, but never “‘certain’’. The main point is that 
our geometric intuition might lead us astray. We could hope to overcome 
this by giving a rigorous definition of what we mean by a surface, and in 
particular a “smooth” or ‘“‘well-behaved” surface. But would such an 
approach be helpful if we wished to extend our methods to functions of 
more than two variables? 


The solution to our difficulties is in fact to abandon the geometric notions 
entirely, and to base our ideas on a few non-geometric axioms. This will 
be one of the tasks of the analysis section of a later course. We shall, as 
it were, use the “‘geometric scaffolding” only as a guide when the building 
starts in earnest. 


Don’t be misled into thinking that the geometric reasoning has been a 
waste of time, for now we have a very clear idea of the results which we 
would like to prove formally. In the future we may well allow our geo- 
metric reasoning to lead the way, but we hope that it will be closely 
followed by rigorous analysis. 


We would like to leave you with one or two thoughts. Could we approxi- 
mate functions of two (or more) variables using an extension of Taylor’s 
series for functions of one variable? Is there any way of deciding if a 
turning point of a function of two variables is a saddle-point? In the 
problem of the aqueduct, is there some other way of shaping the metal 
plates which would give an even greater cross-sectional area? 
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15.4.2 


Discussion 
* 


Unit No. Title of Text 


1 Functions 
2 Errors and Accuracy 
3 Operations and Morphisms 
“ Finite Differences 
5 NO TEXT 3 
Ze Inequalities 
7 Sequences and Limits I 
8 Computing I 
o Integration I 
10 NO TEXT 
ist Logic I — Boolean Algebra 
12 Differentiation I 
13 Integration II 
14 Sequences and Limits II 
15 Differentiation II 
16 Probability and Statistics I 
7 Logic II — Proof 
18 Probability and Statistics II 
19 Relations 
20 Computing II 
21 Probability and Statistics III 
a Linear Algebra I 
23 Linear Algebra II 
24 Differential Equations I 


25 NO TEXT 


26 Linear Algebra III 

es, Complex Numbers I 

28 Linear Algebra IV 

29 Complex Numbers II 

30 Groups I 

31 Differential Equations II 
32 NO TEXT 

33 Groups II 

34 Number Systems 

a5 Topology 

36 Mathematical Structures 
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